Length Distribution of Periodic Orbits 
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The discrete unitary (reversible) analogues of the continuous (irreversible) tent maps are numer- 
ically investigated, in particular, the lengths probability distribution of their periodic orbits. It is 
found that its density can be well approximated by the inverse proportional law. 
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I. INTRODUCTION 

As known, many key features of continuous-time 
chaotic dynamics can be treated in terms of two- 
dimensional or even one-dimensional discrete-time maps 
(see e.g. 0). Among them, most simple and perfect ones 
are the famous Bernoulli map, X.^^i = 2X„modl , and 
the tent map, Xn+i = 1 — |1 — 2X„| , both defined at the 
interval < X < 1. Almost any trajectory of both maps 
is chaotic, that is dense in all this interval. 

In computer, however, any trajectory finishes at zero 
after < D steps only, where D is number of bits under 
use. This contraction of the discretized phase space, from 
total N = 2^ points to single point X = 0, happens 
because the maps are non-invertible, i.e. relation between 
Xn and Xn+i is not an one-to-one correspondence. But 
a suitable slight distortion can make a discretized map 
invertible. Evidently, any such map is nothing but a 
permutation of TV = 2^ states (X„ G [0 : - 1]/A^). 

In particular, the invertible discrete versions of the 
Bernoulli map describe the Mersenne digital auto- 
generators of pseudo-random 0-1-sequences (see e.g. 
For example, the Mersenne map in Fig.l represents tran- 
sitions between 16 states of a particular "toy" 4-bit 
Mersenne generator {D = 4). The general formula of 
such maps is 



Y = 2X + 1- cr{X) , X < N/2 

Y = 2X - N + a{2X - N) , X > N/2 



(1) 



where Y = NXn+i , X = NXn , and any of a{k) , k = 
0,l,..,A^/2 — 1, equals to either zero or unit, a{k) — 
0, 1 (hence this is class of 2^/^ different maps). To see 
corresponding unitary evolution matrix, it is necessary 
to turn the map picture headfirst and insert zeros and 
ones in place of the empty cells and crosses, respectively. 

The map in Fig.l has two immovable points {X = 
and X = 15) and two cycles (periodic orbits) with equal 
lengths L = 7. In applications of much greater Mersenne 
maps (with ~ 70 ^]), the main issue is construction 
of those having most long cycles, with L ^ N . 
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FIG. 1: Example of the discrete Bernoulli map (the Mersenne 
map, see body text). 



Here, we will be interested in modest problem about 
statistics of cycles of discrete maps, to be concrete, defi- 
nite class of discrete tent maps. 



II. CYCLES OF CONTINUOUS TENT MAPS 

First consider the variety of generally asymmetric con- 
tinuous tent maps (CTM) described by 
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X„+i = (l-X„)/(l-a) , X>a, 

where < a < 1. The inversion of (0) reads as 

Xn-i = ^« + (a - Sn)X,, , 5„ = or 1 , (3) 

with S'n = and Sn = ^ representing two variants of 
the inversion. Let us write, with the help of Xq as 
a function of Xl and then equate Xq to Xl at various 
possible values Si,..,Sl = 0,1. Thus we obtain 2^ 
linear equations which determine 2^ points belonging to 
all possible periodic orbits (cycles) of length L . That 
are either irreducible (indivisible) orbits or, if L has a 
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FIG. 2: Two examples of discrete tent maps at 



divider I , composite orbits consisting of L/l repetitions 
of an irreducible orbit with length / . Consequently, 

E = 2^ ' (4) 

1<1 <L 
L mod I — 

where Ni is number of different irreducible periodic or- 
bits with length I , and the sum is taken over all the 
dividers of L . 

In particular, if L is a prime number then I^J yields 

Ni = (2^ - 2)/L {L is prime) , (5) 

since Ni — 2 (there are two immovable points). In gen- 
eral, it is not hard to derive from Q that 

2^ -2(2^^/^^ -[L/2]^ < LNl <2^ -2 (6) 

Here [L/2] = L/2 if L is even, and [L/2] = [L - l)/2 if 
L is odd. At any L » 1, according to lO, 

Nl ~ 2^/L (7) 

At relatively small L, from one finds: A'2 = 1, — 
2, N4 = 3, TVs = 6, Ne = 9, N7 = 18, iVg = 30, 
Ng = 56, Nio = 99, TVii = 186, ... 

III. DISCRETE TENT MAPS 

Next, consider invertible discrete tent maps (DTM) 
representing discrete analogues of the continuous tent 
maps To be concrete, let us introduce the class of 
discrete maps defined by 

y = [m >^<^; 
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D = 4 and a = 1/2 (left) and a = 1/3 (right). 



Here, like in (QJ, X ^ NXn and Y = NXn+i are in- 
tegers, < X,Y < N — 1 ; A is analogue of a in 
< A < N — 1; [x] is an integer closest to x from below 
(with [x] = X if 2; is integer), and {x} means closest 
integer greater than (or equal to) x (that is {x} = x if 
X is integer) . Two examples of maps defined by formula 
(IHl), at = 4, are shown in Fig. 2. 

Let us prove that © prescribers invertible maps, i.e. 
establishes an one-to-one correspondence between X and 
Y. It is sufficient to prove that the upper r.h.s. and lower 
r.h.s. in |(SJ) can not produce one and the same number. 
Indeed, in such a case we would have 

[NX' /A] = Y ^ {N{N - X")/{N - ^)} - 1 , 

with some X' < A and X" > A. This would mean that 

NX' = AY + u, 0<u<A, 
N{N - X") = (TV - A)Y + v , <v < N - A , 

{u and V are some integers), therefore, 

N + X' -X" = Y+{u + v) /N 
< u + v < N 

But, obviously, the latter inequality contradicts the pre- 
ceding equality. The proof is finished. 

It is easy to see also, that absolute value of a devia- 
tion of any DTM © from its CTM prototype lO, with 
a — A/N, (as well as deviation of Q from Bernoulli 
map) does not exceed 1 lowest bit only (±1/A^, in terms 
of Xn). In this sense, any DTM tends to corresponding 
CTM, when — *■ 00. However, the inverted map re- 
mains strongly discontinuous, which is a payment for its 
unambiguity. 

Hence, in essence, the limit of DTM at A^ — > 00 does 
not coincide with corresponding CTM. 
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FIG. 3: (On the left) The length, 
parameter of the map, A , at N = 



Lmax , of most loHg periodic orbit of the discrete tent map, as a function of the asymmetry 
2048 (D — 11). (On the right) The corresponding histogram of Lmax/N. 



IV. CYCLES OF DISCRETE TENT MAPS 

Naturally, there is no simple rule for the cycles (pe- 
riodic orbits) of the DTMs. A number of various irre- 
ducible cycles and lengths of these cycles are extremely 
irregular functions of A and N. In contrary to Q), now 



E 

1<L<W 



LNt=N = 2 



D 



(9) 



where is a number of different irreducible periodic 
orbits with length L (possibly, = 0). For example, at 
D = 12 and some three next values of A this expansion 
looks as 



1 

1*2 
14 



+ 4095 = N , 
+ 2 * 2047 = N 
13*315 = iV, 



where the second multiplier (if any) in each term is num- 
ber of different cycles of a length represented by the first 
multiplier. In other example, for D = 15, 

1*2 + 2 + 4*3 + 8* 30 16 * 2032 = N {A = N/2) , 

1 + 32767 = N {A = N/2 - 1) , 
1*2 + 2 + 4 + 8*5 + 16 * 65 + 24 *2 + 32 + 48 * 9+ 
+60 * 2 + 64 + 72 + 103 + 112 * 20 + 120 + 128 * 10+ 
+ 144 + 176 + 192 * 6 + 224 * 4 + 240 + 1200+ 
+1570 + 8792 + 12999 ^N {A^N/2 + 1) 

Typically, a DTM have a long cycle whose length is 
comparable with the total number of points, A^, i.e. max- 
imal possible length. This is illustrated by Fig. 3 relating 
to D — 11. Quite similar pictures take place also at 
greater D. We see that practically any DTM has a cycle 
with length L > 7V/4 + A^/3. Such long cycles, whose 



length L is comparable with N, can be treated as dis- 
crete analogues of the chaotic trajectories of CTM. 

The left plot in Fig. 3 demonstrates highly irregular 
dependence of the maximal cycle length, Lmaxi on the 
asymmetry parameter A. At the same time, the right- 
hand plot there shows that even slight smoothing of the 
data produces rather regular results. Therefore, it is rea- 
sonable to describe the cycles of DTMs in statistical lan- 
guage, considering some specific subclasses of DTMs in- 
stead of individual maps. 



V. STATISTICS OF CYCLES OF SYMMETRIC 
DISCRETE TENT MAPS 

Let us consider the family of nearly (asymptotically) 
symmetric DTMs determined by the conditions 



N -M 
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N + M 



N,M 



M 
TV 



(10) 



2 - 2 
Since M — > cx), we obtain the growing statistical en- 
semble, which is all the more representative one because 
all the maps defined by H10() do tend to the same limit 
A/N ^ a = 1/2. 

It would be interesting to investigate such the family of 
asymptotically symmetric discrete tent maps (ASDTM), 
in comparison with the usual symmetric CTM. Here, sim- 
plest statistical characteristics of cycles (periodic orbits) 
of the ASDTM will be under our attention. 

Let W{L) designates a density of probability distribu- 
tion of the cycle lengths in this family of maps. In other 
words, practically. 



N 



W{L)^ Nl/Y^Ni 



(11) 
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FIG. 4: (Left plot) The probability, P{L), that an arbitrarily chosen point belongs to some periodic orbit (of the nearly 
symmetric discrete tent maps IIUH ) with lengths less than or equal to L, as a function of L/N, at D = 14 and M ~ 200. 
(Right plot) The binary logarithm of the probability, W{L) , that arbitrarily chosen periodic orbit (of a similar maps family) 
has the length L, via binary logarithm of L, at D = 15 and M ~ 500, in comparison with the analytical estimate of the 
W{L) dependence (the straight line). The inset demonstrates the same comparison for small lengths. 



where is now summary number of periodic orbits 
of length L in aU the maps of the ASDTM family. It is 
useful to introduce also the quantity 



nL)^j^j:iN,, (12) 
1=1 

that is relative (probability) measure of points which be- 
long to all periodic orbits with lengths < L at all the 
maps. 

In reality, with the help of an ordinal PC only, it would 
take rather long time to obtain all the periodic orbits 
if £> > 16 ^ 17. But our computations, performed at 
D = 12^ 16, showed that already D = 13 14 are satis- 
factory values, because next D's increases do not change 
the picture qualitatively (although, of course, providing 
better numeric accuracy). 

Therefore, it is reasonable to prefer calculations at not 
high D {D = 14-^15) but apply slight smoothing over L's 
values. Concretely, the plot W{L) in Fig. 3 represents the 
result of averaging of the exact W{L) (defined by (C3) 
over the intervals [i, L + min((5L, 1)] , with 5 ~ 0.01 . 

Such smoothed probability density, W{L)^ is repre- 
sented by the curved line at right-hand side of Fig. 3. It 
is easy to find that the best fitting for it is nothing but 
the inverse proportional law: 

WiL) . ^ , (13) 



where factor (In N) ^ ensures the normalization, 
Viy(i) « / W{L)dL = 1 

The dependence (|13ll is shown by the straight line, and 
by smooth curve in the inset which demonstrates good 
quality of this fitting for short cycles too. 

The curved line at left side of Fig. 3 shows an example 
of the probability (|12|l . It is not smoothed, therefore, 
formed by many steps with very different heights (like 
the famous devil's staircase). Its closeness to the thin- 
ner straight line, which corresponds to P = i/iV, says 
that the total MN points are distributed approximately 
equally between cycles with different lengths. This is just 
about what the approximation (|13|l says. 

If combining H11() . H12|I . (|13() and the approximation 
P = L/N, we obtain the estimate of mean number of 
cycles per one map of the family: 

1 ^ 

— ^TVi^lniV, (14) 

L=l 

(let us recollect that is summary value for all the AS- 
DTM). At the same time, as the above examples of map 
expansions into cycles do show, fluctuations in number 
of cycles from one particular map to another are sig- 
nificantly greater than the mean value 11411 . According 
to these examples, as well to the W{L) plots, especial 
contribution to the fluctuations comes from cycles whose 
lengths are powers of two. 
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Nevertheless, when rising D from 12 to 16 a definite 
decrease of the fluctuations was noticed. This observa- 
tion pushes us to the hypothesis that formula ifT^ repre- 
sents a true asymptotics of the cycles length distribution 
(fTTjl in the hmit ifTUIl . 

VI. DISCUSSION AND RESUME 

There is a simple naive explanation of the hypothetical 
asymptotics (|13ll . In above derivation of the estimate 
for the cycles numbers of continuous tent maps (CTM), 
the factor 2^ (with L being cycle lengths) arises from 
the two- valued property of their inverse maps (i.e. due to 
their irreversibility) . Since the discrete tent maps (DTM) 
under consideration have univalent inversions, in their 
case this factor must disappear. Thus one deduces the 
inverse proportional dependence of number of the cycles 
(periodic orbits) on their lengths, i.e. comes to H13|l (in 
other words, cycles of different lengths involve approx- 
imately equal amounts of points of the discrete phase 
space). 

In view of such reasonings, the inverse proportional law 
can be expected in case of any sufficiently reach family 
of unitary DTM, not only the ASDTM family defined by 
(|10|l (moreover, in unitary discrete analogues of multi- 



modal piecewise linear CTM, not unimodular tent maps 
only) . 

In factual finite-precision computer simulations of 
more or less general CTM, distinguished from the "patho- 
logical" map Xn+i — 1 — 11 — 2Ar„| (about it see Introduc- 
tion), e.g. asymmetrical CTM Q or the map Xn^i = 
1 - c|l - 2X„| with e.g. c = 0.9999999999.. , it can be 
expected that the law H13|) must manifest itself sooner 
than |7J. Indeed, if the latter rule was true, then the 
very short orbits with lengths L < D only (recollect 
that D is number of bits under operations) would take 
all N = 2^ points of the discrete phase space. In reality, 
with no doubts, almost any choice of an initial point re- 
sults in a long orbit ("chaotic trajectory"), whose length 
L is comparable with iV, while it is hard to hit casually 
into a short orbit. 

To resume, we performed computer statistical analy- 
sis of periodic orbits of unitary (reversible) discrete tent 
maps, and found that probability distribution of the or- 
bits lengths well obeys the inverse proportional law 

I am grateful to Dr. I.Krasnyuk for useful conversa- 
tions. 
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